Compression Induced Folding of a Sheet: An Integrable System 
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The apparently intractable shape of a fold in a compressed elastic film lying on a fluid substrate 
is found to have an exact solution. Such systems buckle at a nonzero wavevector set by the bending 
stiffness of the film and the weight of the substrate fluid. Our solution describes the entire progression 
from a weakly displaced sinusoidal buckling to a single large fold that contacts itself. The pressure 
decrease is exactly quadratic in the lateral displacement. We identify a complex wavevector whose 
magnitude remains invariant with compression. 

PACS numbers: 46.32.+X 46.70.-p 68.60.Bs 81.16.Rf 



Composite structures, containing a fluid substrate cov- 
ered by a thin rigid layer, are commonly found in biolog- 
ical tissues and synthetic coatings. Unlike a freely sus- 
pended sheet, a supported layer has an intrinsic length 
scale arising from the competition of bending and sub- 
strate energy. Thus, e.g., a compressed sheet floating 
on a fluid buckles at a wavelength A = 2'K[B/{pg)]-'-/'^, 
B being the bending stiffness, p the fluid mass density 
and g the gravitational acceleration [l|-l8|. An analogous 
argument holds for an elastic foundation Q. 

In the elastic case, it has long been recognized that this 
extended periodic wrinkling is always unstable against 
localized folding for a sufficiently large system [lO|-ll4|. 
With a fluid substrate the same instability obtains |15l - 
|l9|. Such fold localization has been observed in di- 
verse fluid-supported films — from monolayers and tri- 
layers of nanometer-sized gold particles 15 . 16] , through 
submicron-thick polymer films [ITf, to 10-/zm-thick plas- 
tic sheets 15]. It has been suggested that the localized 
folds, observed in certain surfactant monolayers at the 
water-air interface upon sufficiently fast compression (al- 
beit apparently without prior wrinkling) j2Cll - l25l | , an d be- 
lieved to be important for the function of lun gs |26| . may 
be a manifestation of the same phenomenon |l5l . |l8| . 

The shape of the fold beyond infinitesimal amplitude 
has only been known numerically jlSl . Ilq . The numer- 
ical studies showed puzzling regularities. For example, 
the surface pressure appeared to vary exactly quadrati- 
cally with the displacement. Here we account for these 
regularities by solving the nonlinear equation for the fold 
shape exactly. This allows a much deeper analysis of the 
phenomenon, including its large-deformation limit and 
the point of self-contact, which are central to the fold- 
ing observed in experiments. In a broader context, the 
current work adds an item to the precious collection of 
exactly solvable nonlinear physical problems. 

Consider a thin incompressible elastic sheet of length 
i, width W , and bending modulus B. The sheet is uni- 



axially compressed along the x direction and assumed to 
deform in the xz plane while remaining uniform along 
the y direction; see Fig. [1] Due to incompressibility, the 
configuration of the sheet is completely defined by the 
profile of the angle, 4){s), that the local tangent to the 
sheet makes with the x axis at arclength s. Alternatively, 
we can define a height profile, h{s), where h = sin0 (the 
dot denoting an s-derivative). The region z < /i is occu- 
pied by a fiuid of mass density p = K/g. We focus here 
on localized deformations, and therefore let L -^ oo and 
set 4) — 4> — h — {^dXs^ ±00. 




FIG. 1: Schematic view of the system and its parametrization. 

The energy E = E\^ + E^ contains contributions from 
bending, E^ — (WB/2) J^^ds(p'^, and from the sub- 
strate energy, Es = {WK/2) /^ dsh'^ cos 4>. The dis- 
placement along the direction of compression is 



ds{l 



(1) 



and is related to the pressure by P = dE/dA. For brevity 
we use hereafter units where B = K ~ 1, i.e., we rescale 
energy by B and length by {B/Ky^^, and also \etW = 1. 
The pressure P is scaled by {BKY'^. 

Invoking a dynamical analogy where s stands for time, 
we look for the stable configuration by minimizing the 



action S = J_^ dsC{(j), h, (f>, h), 



£ = -tp'^ + -h'^ costj)- P{l-cosi 



-Q{s){sm^-h), (2) 



where P and Q{s) are Lagrange multipliers replacing, re- 
spectively, the global constraint on A [Eq. ([ij] and the 
local one on the relation between h and (f). (In the case of 
an elastic foundation, the hydrostatic {h'^ /2) cos cj) term 
is replaced by h'^/2 91.) We identify the conjugate mo- 
menta as p0 = dC/ dip = <t> and ph ~ dC/ dh = Q, and 
use them to obtain the Hamiltonian, % = p^ip+phh — C. 
Since C has no explicit dependence on s (the sheet is 
translation- invariant), H is a constant of motion. 



^ = 9P0 + Ph sin (f) - -h^ cos + P(l 



'0)-O, (3) 



where the last equality follows from the boundary con- 
ditions at s — >■ ±00. Equation ([3]) has the consequence 
that, wherever the sheet is horizontal (0 = 0), we have 
\P(f>\ = 1^1, which leads to the geometrical constraint. 



= 0: \<P\ = \h\ ^ \h\. 



(4) 



Hamilton's equation, p^ = —dT-L/dip, yields the follow- 
ing equation of motion: 



' + {h'^/2 + P) sin (j) + Ph cos </> = 0. 



(5) 



Eliminating ph from Eqs. ([3]) and ([5]) and differentiating 
the resulting equation with respect to s, we get 



-H (072 + P)0 + /i = 0. 



(6) 



Equation (jS]) coincides with Euler's elastica problem 
[9|, I27I, [28 1 . It expresses the balance of normal forces 
on an infinitesimal section of the sheet. The last term, 
which usually corresponds to an external normal force 
27j . arises here from hydrostatic pressure. Another dif- 
ferentiation yields the equation in terms of 4> alone. 



"(j)' + [(3/2)02 _^ p]^ + sin0 = 0. 



(7) 



Continuing to assume that (f> and its derivatives vanish 
at infinity, we integrate Eq. ([7]) once to get 



-P02 



1 



i<t>^0. 



(8) 



At first glance, the nonlinear Eq. ([8|) does not seem 
likely to lend itself to a closed-form solution. Inspection 
of the equations above reveals, on the other hand, that 
they readily yield the profile and all of its derivatives at 
s = 0. Let us specialize, for instance, to a symmetric 
deformation about a downward-pointing fold at the ori- 
gin [as in Fig. [2][a)], where 0(0) — 0(0) = 0. We can 
then apply Eqs. ([6|) and ([8]) at s = and, thanks to the 
geometrical condition of Eq. (U), solve for 0(0) and 0(0), 



0(O) = -/i(O) = 2(2-P)i/2, 



and 0(0) = -2(3 - P)(2 - P)!/^. Higher derivatives are 
obtained from Eq. ([7]) and its successive differentiation. 
The complete knowledge of the power-series at s = 
hints that the problem may be integrable. 

Another indication is suggested by the integrable 
physical-pendulum (PP) equation, (p+q'^ sin = 0, whose 
solutions are 



(s) = ±4tan-i(Ae±*«'*), 



(10) 



for any q and A. It is straightforward to show, by in- 
tegrating the PP equation once and differentiating it 
twice, that any of its solutions also solves Eq. (O for 
P = q^ + q~^ + c, where c is an integration constant (set 
hereafter to zero). Thus, Eq. ([T0|) gives complex solutions 
to Eq. (|S]), with the specific complex wavevectors, 

(7 = ±fc±zK, fc = i(2 + P)i/2, ^ = l(2-P)i/2^ (u) 
and an arbitrary amplitude A (the latter following from 



translation invariance). We note that k 



= 1, inde- 



(9) 



pendent of P, while k"^ — k^ = P/2. These solutions 
resemble the 'kink' solutions of the sine-Gordon (SG) 
equation [29|, |30|, albeit in the complex plane. When 
linearized, they coincide with the 'evanescent wave' pro- 
file, which can be stabilized adjacent to a boundary for 
P close to the critical pressure Pc [18[. 

These findings indicate that Eq. ([7]) might bel ong 
to a hierarchy of integrable nonlinear equations [30l . 
3l| . in which the PP equation is a lower-order mem- 
ber. A known hierarchy of equations, referred to as 
the stationary-sine-Gordon- modified-Korteweg-de-Vries 
hierarchy [3l[, indeed contains the stationary SG equa- 
tion, the PP equation, and Eq. ([7|) as the first, sec- 
ond, and third members, respectively. To our knowledge, 
equations in this hierarchy beyond the PP equation have 
not been linked before to physical phenomena. 

To construct localized real solutions out of the complex 
ones given in Eq. ([rU]) . we borrow a scheme from the SG 
problem. In 'light-cone' coordinates [u = (x + t)/2, v = 
{x — t)/2], the SG equation, dxxfp ~ dtt(p ~ duv4> — sin0, 
is invariant to the scaling u -^ qu, v — > v/q by an ar- 
bitrary scale factor q. Given three known solutions of 
this equation, 0j (j = 0,1,2), one can construct an- 
other solution, 03, using the imphcit 'ladder' rule [30|, 
tan[(03 - 0o)/4] = [(gi -I- q2)/{qi - 92)] tan[(0i - 02)/4], 
where gi and (72 are arbitrary scale factors for 0i and 
02. We attempt the same procedure while substituting 
for qi and 52 two of the specific P-dependent wavevec- 
tors found in Eq. pT]) . Choosing 0o = 0, 0i = 
4tan^i(e'«i"), and 02 = Atan-^ie''^'''), with qi = k - in 
and 52 ~ —k—in, we obtain the odd function tan(03/4) = 
(k/A;) sin(fcs)/cosh(Ks). Using instead tan(0i/4) = ie"^^" 
and tan(02/4) = — ie*''^'*, we get the even counterpart. 
Substitution of these two functions in Eq. ^ confirms 
that they indeed solve it exactly. Thus, the following are 



exact localized shapes of the angular profile: 



symmetric fold: (j>{s) = 4 tan 
antisymmetric fold: (f>{s) = 4tan~ 



K, sin(fcs) 

fccosh(Ks) 

KCOs{ks) 



kcosh{K,s) 



(12) 



These functions match the 'breather' solutions of the SG 
equation [29[ when those are projected onto the light 
cone, {s — X = t = u,v — 0). Due to the symme- 
tries under reflection about the z axis, reflection about 
the s (or x) axis, and translation along s, the functions 
±(/)(±s + So), where (/)(s) is either one of the functions 
in Eq. (J12p and sq an arbitrary constant, are solutions 
as well. The existence of odd and even solutions then 
follows from the aforementioned ladder rule. 

Evidently, the equations simplify when n becomes 
small, as P ^ Pc = 2; then, e.g., the symmetric fold has 
(/> ~ 4Ksins/ cosh(Ks), which itself is vanishingly small. 
This is the regime of incipient buckling discussed previ- 
ously [l|, y, la UM, UM ■ The buckling is always localized, 
but the localization length diverges as the threshold is 
approached |18l.ll9j. 

The solution implies very simple relations among the 
pressure P, the displacement A, the central height |/i(0)|, 
and the energies. The decay parameter is exactly linear in 
the displacement, k = A/8. Indeed, the expressions for 
k and k [Eq. (jlip ] exactly match the complex wavevector 
obtained from a linear analysis of the 'evanescent-wave' 
for P — > Pc [18|. Consequently, the pressure is exactly 
quadratic in A: P — 2 — A^/16, as previously deduced in 
that limit ^1&] . The maximum amplitude of a symmetric 
deformation is \h{0)\ = A/2. The bending and substrate 
contributions to the energy are E\^ = A and £"3 = A(l — 
A^/48). We note that the energies and pressures are 
identical for the symmetric and antisymmetric cases. 

Figure [2] shows the progression of the symmetric 
and antisymmetric folds as the lateral displacement in- 
creases and the pressure decreases. The configurations 
have been calculated from Eq. (IT^ according to the 
parametrization: x{s) = L ds' cos (f>{s'), z(s) = h{0) + 
Jjj ds' sin (j>{s'), where h{0) is given for the symmetric fold 
by Eq. (J9]) and for the antisymmetric one by h{0) = 0. 

The symmetric fold is found to contact itself at a 
small positive pressure, P ~ 0.040 (corresponding to 
A ~ 5.6 ~ 0.89A). Self-contact of the antisymmetric fold, 
by contrast, requires a substantial negative pressure (i.e., 
tension) of P 2± -0.70 (A ~ 6.6 ~ 1.05A). Thus, in the 
case of an antisymmetric configuration the stress in the 
sheet vanishes prior to self-contact. One can examine the 
solutions also beyond self-contact, where they produce 
self-intersecting configurations which are unphysical for 
a sheet. In particular, at P = —2 the oscillations in cj) 
disappear. These configurations are shown in Fig. |3l 

These remarkably simple yet exact results are conse- 
quences of the high level of symmetry characteristic of 



integrable nonlinear problems [30| . The case of an elas- 
tic foundation [10| is physically less simple, since in-plane 
shear forces (and not merely normal hydrostatic ones) are 
exerted on the sheet. Although this problem is found to 
obey the same constraint at extrema [Eq. (JH)], it does 
not exhibit the regularities described above. 

The solution presented here provides precise knowl- 
edge of shapes and energies in a large class of wrinkling 
and folding systems. It enables a new means for mak- 
ing precise force actuators and transducers on any scale 
where uniform, thin sheets can be made, for example. 
The solution improves the prospects for understanding 
the unstable motion resulting from folding |23l425| and 
the observed buckling of nanoparticle monolayers into 
trilayers [l6|. At the molecular scale, it provides a start- 
ing point for quantifying the effects of compressibility and 
self-attraction of surfactant monolayers, as well as the in- 
fiuence of non-fluid aspects of the substrate. More basi- 
cally, compression-induced folding appears to be a previ- 
ously unrecognized integrable solitary wave phenomenon, 
like the sine-Gordon chain and the Korteweg-de Vries 
hydrodynamic soliton. The results above may be used 
to construct more complex, multiple-fold shapes. The 
fundamental reason for the integrability of the problem 
remains to be understood. It is to be hoped that this 
understanding will reveal a broader class of physical sys- 
tems which are integrable for the same reason. 
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